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The energy band structure and energy splitting due to quantum tunnel- 
ing in two weakly linked Bose-Einstein condensates were calculated by using 
the instanton method. The intrinsic coherent properties of Bose Josephson 
junction were investigated in terms of energy splitting. For Ec/Ej <C 1, the 
energy splitting is small and the system is globally phase coherent. In the 
opposite limit, Ec/Ej 3> 1, the energy splitting is large and the system be- 
comes a phase dissipation. Our reslults suggest that one should investigate 
the coherence phenomna of B J J in proper condition such as Eq / Ej ~ 1 . 
PACS numbers: 03.75.Fi, 05.30.Jp, 64.60.My, 67.40.Fd 

I. INTRODUCTION 

Two weakly linked Bose-Einstein condensates (BECs) behave as the superconductor 
Josephson junction. This remarkable feature has been investigated both theoretically [1-6] 
and experimentally [8-12]. The existence of a Josephson current through a potential barrier 
between two superconductors or between two superfmids is a direct manifestation of macro- 
scopic quantum coherence [13]. The experimental realization of Bose-Einstein condensation 
(BEC) of weakly interacting alkali atoms [8] has provided a route to study neutral super- 
fluids in a controlled and tunable environment [9]. The possibility of loading a BEC in a 
one-dimensional periodic potential has allowed the observation of quantum phase effects on 
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a macroscopic scale such as quantum interference [10], superfluidity on a local scale [11] and 
an oscillating atomic current in Josephson junction arrays [12]. 

The analogy of the volt age- current characteristic in superconductor Josephson junction 
was proposed theoretically in Bose- Josephson junction [2-5]. The macroscopic BEC's co- 
herence has been demonstrated by interference experiments [14], and the first evidence of 
the coherent tunneling in an atomic array, related to the "ac" Josephson effect, has been 
reported [10]. A "dc" current-biased Bose- Josephson junction can be simulated by a tun- 
neling barrier moving with constant velocity across the trap [3] . At a critical velocity of the 
barrier (proportional to the critical tunneling current), a sharp transition between the "dc" 
and "ac" Josephson regimes was predicted. Thus two weakly linked condensates exhibit 
the analog of the resistively shunted superconductor Josephson junctions. The "secondary" 
quantum macroscopic effects in small capacitance Josephson junction had attracted much 
attention both from theorists and experimentalists for decades [15-18]. These phenomena 
manifest the quantum behavior of a Josephson junction as a macroscopic object, in contrast 
with such "primary" quantum macroscopic phenomena as the Josephson junction itself [15]. 
In these cases, quantum fluctuations of the phase difference ip across the junctions become 
important. This necessitates treating the phase as a quantum operator ip, which is canoni- 
cally conjugate to the operator N. Novel macroscopic quantum phenomena, such as "Bloch" 
oscillaition, had been reported in current-biased Josephson junctions [15]. To understand 
these effects and present a general picture of the low temperature dynamics of Josephson 
junctions, a simple theory [16-18] has been suggested based on the extended coordinates, 
if E [— oo, +oo]. 

More recently, the first experimental of BJJ have been realized in a purely magnetic 
double- well potential [24]. It will stimulate the further study on the BJJ. Then, can we 
investigate the "secondary" quantum macroscopic effects [15] in two weakly linked BECs? 
It was argued in [2,3,6] that the dynamics of the phase of this system can be mapped onto the 
sine-Gordon quantum mechanical Hamiltonian (1). In this paper, we calculated the energy 
band structure, "Bolch" wave and the energy splitting of the Bose- Josephson junction due 
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to the quantum tunneling using the instanton method [21]. As a simply application, we 
investigate the coherent properties of the Bose-Josephson junction. The results from this 
method agree exactly with the prediction of Stringari [6]. 

II. MODEL AND ITS INSTANTONS SOLUTION 

All quantities describing the junction should be considered as operators rather than the 
classical variables. The operators corresponding to the main variables, the phase difference f 
and the particle number N of the Bose-Josephson junction, satisfy the commutation relation: 
[ip, N] = i, so that p and N are related by the Heisenberg uncertainty relation which was 
disregarded by the "classical" theory of the Bose-Josephson junction [2]. 

The idealization of two weakly linked BECs can be described as a two-mode bosonic 
system [2]. In the "phase" representation the relevant quantum observables are the difference 
of phases and number of atoms between the two condensates in each trap. The Hamiltonian 
can be written in terms of a "quantum pendulum" (Mathieu) equation 

The "charging energy" E c and the "Josephson coupling energy" Ej can be calculated as 
overlap integrals [2,3], 

E c = 2gJ dr<f>i(r) = 2g J dr<f> 4 2 (r), 

Ej = Nj rfr$Kr)[-|^V 2 + V ext + ^($?(r) + ^(r))]$ 2 (r) (2) 

with the one-body wave functions $i(r), $ 2 (?") localized in the trap 1, 2; and 
/ dr$l(r)$ 2 (r) = 0, / dr$l 2 (r)$ 1:2 (r) = 1, g = ^ SL \ a is the scattering length and m 
the atomic mass; N = a\ai + is the total number of atoms. 

For convenience, we modify the potential as Uj(<f>) = Ej(l + cos0). Then, the Hamil- 
tonian H in Eq. (1) can be considered as one-dimensional quantum particles with mass 
rricf, = Jrj moving along the axis in sine-Gordon potential Uj(<p). Due to the translational 
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symmetry (0— > + 27r) of the Hamiltonian (1), the set of the eigenfunctions should include 
the Bloch wave functions 



(0) 



u ntq (0) exp (iq<f>) , 



n = 0,l,2,..., 



Un,q (0) 



U ntq (0 + 27T) , 



— oo < q < oo 



(3) 



where q is an arbitrary (real) constant vector. Substitution of the wave function (3) in the 
Schrodinger equation leads immediately to the picture of band energy spectrum and relative 
well known effects in solid state theory. From Ref. [17,18], this is nothing but Mathieu equa- 
tion, so that these functions can be readily calculated. Some of their asymptotic properties 
(in tight-binding limit and weak-binding limit) can be expressed analytically. 

The energy band structure and energy splitting for the sine-Gordon potential can be 
calculated alternatively by instanton method [19,20]. The advantage of this nonperturbative 
method, as presented here, is that it gives not only a more accurate description of the 
tunneling phenomena but also a comprehensive physical understanding in the context of 
quantum field theory. 

The effective Lagrangian is 



The classical solution which extremizes the action is seen to satisfy the equation of motion 



where the Wick rotation t — it has taken the system into Euclidean time. Equation (5) can 
be regarded as the equation of motion of a pseudo-particle with the classical energy E d > 0, 
which is a constant of integration. With E c i being confined to a region < E c i < Ej, the 
configuration C becomes periodic such that c (r + T) = c (r), which now corresponds to 
the periodic boundary condition in the space coordinate [19]. The classical solution is 




(4) 




(5) 



C = 2 arcsin[/csn(u;o(T + To))]) 



(6) 
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where ojq = y/EjE c /h is the classical plasma frequency [6], sn is the Jacobian elliptic function 
with modulus k = yjl — E c i/2Ej. The elliptic function sn(w (' r + To)) has a period 4)C(k) 
with K{k) the complete elliptic integral of first kind. For zero energy E d — (k — > 1), the 
periodic solution reduces to the vacuum instanton configuration <p c — > 2 arcsin[tanh(wo(T + 
to))]- 

III. ENERGY BAND STRUCTURE AND THE TRANSITION AMPLITUDE FOR 

QUANTUM TUNNELING 

The sine-Gordon potential has an infinite number of degenerate vacua. Quantum tunnel- 
ing between neighboring vacua leads to the level splitting, while the levels extend to bands 
due to the translational symmetry expressed by Uj(<j) + 2n) = Uj((f>). In the narrow-band 
approximation one finds, for the energy, the expression 

E = Ei + ]T J(R m - R n )e^ R ™- R "\ (7) 

n 

where E^ denotes the ith eigenvalue of the energy in each well for the harmonic oscillator 
potential U(<f>) = |(0 — R n ) 2 , R n = 2n7r is the position of the nth minimum, J(R n r — 
R n )= J ip* (<£>— R n ')[Uj((j)) — U (0)]^(0— R n )d(f) is the overlap integral, ipi is the eigenfunction 
corresponding to eigenvalue E iy and q is the Floquet parameter associated with the Bloch 
wave function. If only the contribution from the nearest neighbors is taken into account, 
i.e., J(R n i — R n ) for | n' — n\ > 1 is taken to be zero, the energy band formula reduces to 

E = Ei + 2Jcos(27rg). (8) 

The parameter J is just the level splitting resulting from quantum tunneling (the wave 
functions are periodic for q = 0). We will consider the case of potential wells surrounded 
by very high potential barriers with correspondingly small tunneling contribution to the 
eigenvalues. They are almost those of degenerate harmonic oscillators, and in this asymptotic 
case we are not concerned with the entire bands but only with their edges which correspond 
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to alternately even and odd states. Then, we suppose \ i)n, \ i)l are degenerate eigenstates in 
neighboring wells, respectively with the same energy eigenvalue Ei such that H° | i)r,l = I 
i)r,l where H° is the Hamiltonian of the harmonic oscillator as the zero order approximation 
of the system. The degeneracy will be removed by the small tunneling effect which leads to 
the level splitting. The eigenstates of the Hamiltonian H become 

I i)o = ^= (I i)n- I i) L ) , I i)e = ^= (I i)n+ I (9) 

with eigenvalues Ei ± AEi, respectively. AEi denotes the shift of one oscillator level. It is 
obvious that r(i \ H — H° \ i) L = 2 J = AEi. In the following we calculate this energy shift 
AEi as resulting from periodic instantons and instanton-anti-instanton pairs. 

The amplitude for a transition from one well to its neighboring well at the energy Ei due 
to instanton tunneling can be written as 

A + ,- = + (^|e- 2 %|^)_ ~ ±e- 2 ^sinh(2^), (10) 

where we neglect overlap of the wave functions which dominate over either well. The am- 
plitude (10) can also be calculated with the help of the path-integral method, A + _ = 
J V'e,h-(0/)^(0/) r /j < t>ii T i)'>pE-{. < t>i)d<t>fd(i)i 1 where the Feynman kernel is defined as usual by 

K(<j> f ,T f ;<f> i ,T i ) = /%[0] e -^ (11) 

■>4>i 

with (fij = 4>{jf), (pi = 4>{ji) and tj — t-i = 2T. What we are interested in is this expression in 
the limits fa — > —a, <pf a (±a are the turning points), namely, the tunneling propagator 
through one of the barriers. S = f£ [\ m <t>{^f) 2 + Uj{fa)]dT is the Euclidean action of the 
pseudo-particle, and ipE^+ifa) (i>Ei,-(fa)) is the wave function of the right- (left-) hand wells. 

The functional integral K(<f>f,Tf;fa,Ti) can be evaluated with the stationary method 
by expanding (r) about the classical trajectory C (r) and thus we set (r) = C (r) + 
x(t), where x ( r ) is the small fluctuation with boundary conditions x( T i) = xijf) — 0. 
Substitution of <fi (r) into Eq.(ll) and keeping only terms containing x ( r ) U P to the one- 
loop approximation yields K(<j>f,Tf;fa,Ti) = exp[— where I = J*^2o° T^[x\ e ~ SS ^ h 
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is the fluctuation function integral with the fluctuation action SS = f^ f xM\dr, where 
M = —^pr + V"(0 c (t)) is the second variational operator of the action. The classical 
action S c (r) is evaluated along the trajectory C (t) so that 

I d<p c (r) 



2 V dr 



dr 



= 2 ET + —[S(k)-k' 2 JC(k)}, (12) 

where £(k) denotes the complete elliptic integral of the second kind and k' 2 = 1 — k 2 . 
Following the standard procedure of the periodic instanton calculation in Refs. [19] and 
[20], the functional intergal / can be wrtten as / = ^7= N (r;) N (t/) J^ f N d 2 T ^ , where 
N (t) = d<t> ^ ; is the zero eigenmode of M. To obtain the desired results, the contributions 



from the infinite number of the instantons and anti-instanton pairs have to be taken into 
account. The total amplitude is found to be 



sinh( ac(Fr T) ' < i3) 

with W = ^[£{k) - k ,2 K(k)}. Comparing (13) and (10) leads to AE = i ^exp(-f ). 
Rescaling this formula in the unit of Ej, we arrive at our final result. 



f = \/f4^WI (£W -*' 2KW)l ' <14) 
The formula (14) shows that the widths of the energy bands are very sensitive to the di- 
mensionless parameter of E c /Ej, which is a critical parameter characterizing the dynamical 
properties of the system such as the coherence [6] and the transition [15]. A direct appli- 
cation of this energy splitting is to investigate the coherence properties of BJJ. Quantum 
coherence requires that the relative phase of the order parameter should be preserved over 
times of the order which is just the tunneling amplitude [23]. From the above calcu- 
lation, we know that the energy splitting AE is proportional to T. Therefore the energy 
splitting describes the coherent properties of the system. 

Fig. 1 shows the energy splitting AE/Ej as a function of the ratio Ec/Ej in the 
ground state (E c i = 0). The figure shows that for values of Ej smaller than Ec the energy 



splitting is significantly increased, indicating the occurrence of a continuous transition to 
the phase dissipation (the number squeezed regime) [6,7]. In the limit E c /Ej <C 1 the 
system undergoes small oscillations around the equilibrium. In this limit the correlation 
between the neighboring wells is small. One can regard the systems as a globally coherent 
object described by a unique order parameter. On the contrary E c /Ej 3> 1 the behavior of 
the system is very different. The quantum fluctuation is enhanced due to the increasing of 
the tunneling between the neighboring wells in phase space, showing that the relative phase 
between two condensates are distributed in a random way. At the same time, the fluctuation 
of the relative number of atoms in two traps becomes smaller and smaller. 

It is interesting to investigate the properties of the systems under the excited states 
(Ed 7^ 0) which may be excited by the thermal fluctuations or other reasons. We plotted 
the energy splitting as a function of the dimensionless parameter E d /Ej for two different 
values of E c /Ej in Fig. 2. One clearly observes that even if quantum effects are small the 
decoherence due to higher excited states may become important. It is interesting that the 
respone of the system to the decoherence perform a total different behavior. States with 
small quantum effects are more sensitive to the decoherence fluctuations than those with 
relatively large vaules of E c /Ej. The system with little higher E c /Ej state should be a 
good candidate to investigate the quantum coherence phenomena, in Ec/Ej ~ 1. Because 
in this case the system with a good coherence will preserve its coherence for a relative large 
range of the energy perturbation. 

Energy spectrum and energy splitting due to quantum tunneling in BJJ have been cal- 
culated by means of instanton method. Based on this energy splitting formula, we also 
investigated the coherence property of BJJ. Our results agree exactly with that in Ref. 
[6]. This analysis makes it possible to investigate the "secondary" quantum phenomena in 
BJJ (see, for instance, [18]) and presents a general picture of low temperature dynamics of 
Josephson junctions in BEC. 
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Figure Captions 

Fig. 1. Ground state energy splitting as a function of the ratio Ec/Ej. 
Fig. 2. Energy splitting as a function of the ratio E/Ej for E c /Ej = 1 (solid line), 
E c /Ej = 3 (dashed line). 
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Fig.2 



